Structure and optical properties of high light output halide scintillators 
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Structural and optical properties ol several high light output halide scintillators and closely related 
materials are presented based on first principles calculations. The optical properties are based on 
the Engel-Vosko generalized gradient approximation and the recently developed density functional 
of Tran and Blaha. The materials investigated are BaBr2, BalBr, BaCb, BaF2, Bal2, BH3, Cal2, 
Cs 2 LiYCl 6 , CsBa 2 Br 5 , CsBa 2 I 5 , K 2 LaBr 5 , ^LaCls.I^Lals, LaBr 3 , LaCl 3 , SrBr 2 , and YI3. For 
comparison results are presented for the oxide CdWCU. We find that the Tran Blaha functional 
gives greatly improved band gaps and optical properties in this class of materials. Furthermore, we 
find that unlike CdWC>4, most of these halides are highly isotropic from an optical point of view 
even though in many cases the crystal structures and other properties are not. This general result 
is rationalized in terms of halide chemistry. Implications for the development of ceramic halide 
scintillators are discussed. 

PACS numbers: 78.20.Ci,78.20.Bh,61.66.Fn 
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I. INTRODUCTION 

Scintillators are widely used in radiation detection ap- 
plications including medical imaging, oil well drilling, 
nuclear security and high energy physics experiments. 
These materials function by emitting light when excited 
by ionizing radiation, such as gamma rays. This light is 
then coupled to a photomultiplier or other light detector 
to produce electrical signals. The performance of scin- 
tillators is characterized by their light output, normally 
given in terms of photons per MeV of excitation energy, 
proportionality (how linear the light output is as a func- 
tion of excitation energy), density (related to stopping 
power), decay time and energy resolution. 1 

Scintillator performance depends on energy transport 
between the energy absorption events and the scintilla- 
tion centers, suppression of non-radiative recombination 
channels and fast efficient light emission from the scintil- 
lation center. These scintillation centers can be intrinsic, 
via a mechanism such as decay of a self trapped exci- 
ton, or, as is commonly the case, at an activation center 
(e.g. Ce 3+ ions substituting for La 3 + in LaBra) where 
electron hole pairs generated by the radiation recombine. 
These processes are sensitive to the details of the mate- 
rial and in the case of materials with activators the envi- 
ronment of the activator ions. Halides generally have 
soft lattices that disfavor non-radiative recombination 
and can typically incorporate rare earth and other ac- 
tivator ions. More importantly, halide chemistry is very 
rich, with a wide variety of crystal structures that pro- 
vide different bonding topologies for energy transfer and 
various environments for activator ions^i^ In fact, one of 
the best known scintillators is the halide, Nal activated 
with Tl + . This material has high light output, but is 
slow, non-proportional and has poor energy resolution. 4 
The finding that Ce 3+ activated LaBr3 is a very high 
light output proportional scintillator with energy reso- 
lution better than 3% at 662 KeV- has led to renewed 



interest in halides as scintillator hosts, especially for spec- 
troscopic gamma ray detection. This interest has resulted 
in the discovery of several other interesting halide scin- 
tillators, including heavily Eu 2 activated Sri2, which is 
very proportional and has a light output exceeding that 
of LaBr3j£ and Ce 3+ activated YI3, which is another very 
high light output material^ 

Full characterization of the optical properties of scin- 
tillators is very useful both from the point of view of im- 
proving the design of systems as regards light coupling 
and also importantly in selecting candidate materials for 
ceramic scintillators. However, full optical characteri- 
zation of these halides is complicated by sensitivity to 
moisture and other experimental difficulties. As a result 
only limited data is available. 

The key requirements for a ceramic scintillator are sin- 
terability and optical isotropy. This latter requirement 
comes from the need to avoid light scattering at mis- 
oriented grain boundaries. It is commonly thought that 
because of this requirement cubic material is needed for a 
true transparent ceramic. However, radiation detection 
is not an optical imaging application. Therefore weak 
scattering and image distortion are not as detrimental as 
they would be in an optical application such as a ceramic 
lens. In fact ceramic scintillators based on monoclinic 
Lu 2 Si05 (LSO) have been demonstrated^ We showed 
in previous work that the high light output halide scintil- 
lator, Srl 2 is in fact very nearly optically isotropic in spite 
of its orthorhombic (Pbca) crystal structured This was 
an expected result considering the strongly orthorhom- 
bic lattice. It is of interest to determine whether this 
is the case for other high light output halide scintilla- 
tors. Here we present a consistent set of first principles 
data for structural and optical properties of a number of 
high light output halide scintillators and closely related 
materials. We find, quite unexpectedly, that the halides 
we investigate have remarkably little optical anisotropy, 
even though a number of them are very anisotropic from 
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other points of view. 



II. APPROACH 

The density functional calculations reported here were 
performed using the general potential linearized aug- 
mented planewave (LAPW) method as implemented in 
the WIEN2K codeJ^- We used well converged Brillouin 
zone samples and basis sets, with the standard LAPW 
augmentation plus local orbitals, 12 Relativity was in- 
cluded at the scalar relativistic level except for Bil 3 , 
where spin-orbit was included for the electronic struc- 
ture. 

The crystal structure plays a fundamental role in deter- 
mining the electronic and optical properties of a material. 
We began our calculations by fully relaxing all free inter- 
nal atomic positions consistent with the crystal symme- 
try for each material. The lattice parameters were held 
fixed at their experimental values, which are no doubt 
more precise than the values that can be obtained using 
density functional calculations. This relaxation was done 
using the generalized gradient approximation (GGA) of 
Perdew, Burke and Ernzerhof (PBE).—. 

The PBE functional, like other standard generalized 
gradient approximations, is based on the total energy 
in terms of the coupling constant averaged exchange 
correlation hole and is designed to reproduce the total 
energy, 14 While these functionals are very useful in ob- 
taining structures and other properties related to total 
energies, they underestimate, often strongly, the band 
gaps of most semiconductors and simple insulators. Ac- 
cordingly, for the optical properties we use two other 
functionals. The first is the Engel-Vosko GGA (EV). 15 
This functional was designed to reproduce the exchange- 
correlation potential rather than the total energy, and 
gives improved band gaps»i& We used it in our prior 
studies on transport in PbTe and optical properties of 
Srl 2 and Bi4Ge30i2 scintillators! 10 i 17 ' 18 The second is 
the semi- local functional of Tran and Blaha (TB-mBJ). 19 
This relatively recent functional is more sophisticated 
than the Engel-Vosko GGA and has been shown to give 
very much improved band gaps for a variety of semicon- 
ductors and insulators. We find that where comparison 
with experiment is possible this functional also gives very 
much improved band gaps as well as optical properties for 
these halides. As such, we focus on the results obtained 
using the TB-mBJ functional. The electronic structures 
were calculated with these two functionals based on the 
relaxed crystal structures from the PBE calculations. 
Optical properties were then obtained using the dipole 
matrix elements with the WIEN2K optical packaged A 
0.1 eV broadening was applied to the spectra. 



III. MATERIALS AND STRUCTURES 

We begin with the calculated structural parameters 
and brief introductions to the materials that we study. 

BaF2 is a very well characterized material. It has been 
applied as a scintillator, both in pure form and with 
activation<22r— In pure form it has a very fast compo- 
nent (~0.8 ns), as well as a slow component, which can 
be at least partially suppressed by La doping. Fast re- 
sponse is of importance in applications with very high 
count rates or where timing is critical. It occurs in the 
cubic (Fm3m) CaF2 structure, Ba on Wyckoff site 4a 
(0,0,0) and F on 8c (1/4,1/4,1/4). In our calculations we 
used the experimental lattice parameter^ a=6.2001 A. 

Several of the alkaline earth di- halides, AeX2, are high 
light output proportional scintillators when activated 
with Eu 2+ £>2£-22 As with BaF 2 , activation with Ce 3+ 
is also possible in some cases, 30 Undoped BaCl2 has a 
very fast response, with a short component lifetime of 1.6 
ns£L SrI 2 :Eu is equal to or superior to LaBr 3 :Ce 3+ as re- 
gards light output and proportionality, although it suffers 
from a slower response^ Even though the material is or- 
thorhombic, we found it to be optically nearly isotropic J£ 
Cal2:Eu 2+ is another material with very high light out- 
put, ~ 100,000 photons/MeV&2Z BaIBr:Eu 2+ crystals 
have been shown to have a light output of 81,000±3,000 
photons/MeV with a 662 KeV energy resolution bet- 
ter than 5%M Regarding structure, BaCl2, BaBr 2 , Bal2 
and BalBr occur in an orthorhombic Pnma, PbCl2 type 
structure with four formula units per cell. The structural 
parameters as determined from relaxation are given in 
Table HI SrBr2 has a large tetragonal unit cell {PA/nz) 
with ten formula units per celli^ 3 . The calculated inter- 
nal parameters are as given in Table [II] Cal2 is hexag- 
onal (P3ml), with Ca on site la (0,0,0) and I on 2d 
(1/2, 2/3, z). In our optical calculations we used the ex- 
perimental lattice parameters, a=4.49 A, c=6.975 A, 
with the calculated z=0.2535 from total energy mini- 
mization with the PBE GGA (the reported experimental 
value is 0.25)^ 

LaCl3 and LaBr3 activated with Ce 3+ are very high 
light output, proportional scintillators with excellent en- 
ergy resolution^ They are hexagonal (P63/777,) with the 
UCI3 structure type. There are two formula units per 
cell, with La on site 2c (1/3,2/3,1/4) and the halogen 
on Qh (u,t>,l/4). We used the experimental lattice pa- 
rameters, a=7.4779 A, c=4.3745 A, for LaCl 3 (Ref. m) 
and a=7.9648 A, c=4.5119 A, for LaBr 3 (Ref. [H). The 
calculated internal parameters were m=0.9145, i>=0.6132 
for LaCl 3 and u=0.9144, v=0.6159 for LaBr 3 . 

Cs2LiYCl6 is a member of the elpasolite family. The 
elpasolites are cubic Fm3m halides, with general formula 
AI2A2RX6, with Al and A2 alkali metals, X a halogen 
and R a rare earth or other trivalcnt clement. Al is on 
site 8c (1/4,1/4,1/4), A2 is on 46 (1/2,1/2,1/2), R is on 
4a (0,0,0) and X is on 24e (m,0,0). A large number of 
such compounds are known^ However, only a fraction 
have been studied as potential scintillator materials 42r— 
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TABLE I: Structural properties of orthorhombic Prima PbCk 
type, BaCb, BaBr2, Bal2 and BalBr. The lattice parameters 
are from experiment (Refs. H3, 13 and |49T) . while the internal 
parameters are from total energy minimization using the PBE 
functional. For BalBr, XI is I and X2 is Br. All atoms are 
on Wyckoff site 4c with y=0.25. 





BaCl 2 


BaBr 2 


Bal 2 


BalBr 


a (A) 


7.878 


8.276 


8.922 


8.6 


6(A) 


4.731 


4.956 


5.304 


5.12 


c(A) 


9.415 


9.919 


10.695 


10.31 


Ba x 


0.2481 


0.2470 


0.2466 


0.2302 


Ba z 


0.3838 


0.1163 


0.1166 


0.1228 


XI x 


0.1423 


0.6426 


0.6424 


0.9722 


XI z 


0.0708 


0.0714 


0.0721 


0.6701 


X2 x 


0.5274 


0.5266 


0.5246 


0.3538 


X2 z 


0.8313 


0.6686 


0.6665 


0.4325 



TABLE II: Internal structural coordinates of tetragonal 
PA/nz (No. 85) SrBr2 as determined by total energy min- 
imization. The lattice parameters were fixed at their experi- 
mental values of a=11.630 A, and c=7.146 A. (Ref. \M) 





Sri 


Sr2 


Brl 


Br2 


Br3 


Br4 


X 


0.4141 


0.7500 


0.5402 


0.5426 


0.2500 


0.2500 


y 


0.6035 


0.7500 


0.6527 


0.8386 


0.7500 


0.7500 


z 


0.7520 


0.1487 


0.3761 


0.9004 


0.0000 


0.5000 



CS2L1YQ6 is of particular interest because of its high Li 
content, which makes it useful as a neutron detector. We 
used the experimental lattice parameter for Cs 2 LiYCl6 
from Reber and co-workers^ 4 - a=10.4857 A, with the cal- 
culated internal parameter u=0.2517 (an experimental 
value of 0.25046 was reported by Reber and co-workers). 

The scintillation properties of Ce 3+ activated 
K^LaCls, K2LaBr5 and K^Lals were investigated by 
van Loef and co-workersj 45 ' 46 These isostructural or- 
thorhombic compounds showed high light output, up to 
55,000 photons / MeV (for K 2 LaI 5 :Ce 3+ ). The iodide 
also showed a reasonable decay time of 24±5 ns and 
662 KeV energy resolution of 4.5±0.5%, although use 
of this scintillator is complicated because of self-activity 
associated with K. Our calculated structural properties 
are give in Table IIIII 

CsBa2Br5 and CsBa2l5 occur in a monoclinic P2±/c 
(No. 14) structure.—^ 3 - Eu 2+ doped CsBa 2 l5 was re- 
ported to have a light yield of approximately 97,000 pho- 
tons per MeV and to be less hygroscopic than LaBr3 or 
Srl2<££ The calculated structural parameters are as given 
in Table El 

YI3 and B1I3 are related materials. Ce 3+ activated YI3 
was recently reported to have a very high light yield of 
almost 100,000 photons per MeVZ^l The related com- 
pound B1I3 is of interest mainly as a semiconductor for 
radiation detection^ rather than as a scintillator. We 
include it here for comparison with YI3. Unlike YI3, 
Bil 3 has substantial covalency between Bi 6p and halo- 
gen p states, which leads to enhanced Born charges and 
may be expected to affect the optical properties,— These 



TABLE III: Structural properties of orthorhombic Prima 
K^LaAs, A=Cl,Br,I. The lattice parameters are from exper- 
iment (Refs. [Ho and l5ll) , while the internal parameters are 
from total energy minimization using the PBE functional. 
The y coordinate is not listed for the atoms that are on Wyck- 
off site 4c for which y=0.25. 





K 2 LaCl 5 


K 2 LaBr 5 


K 2 LaI 5 


a (A) 


12.742 


13.36 


14.332 


6(A) 


8.868 


9.272 


9.912 


c(A) 


8.022 


8.462 


9.132 


K x 


0.3279 


0.3277 


0.3264 


Ky 


0.9952 


0.5049 


0.5035 


K z 


0.0475 


0.4465 


0.4419 


La x 


0.0066 


0.0070 


0.0064 


La z 


0.0774 


0.4202 


0.4182 


XI x 


0.1818 


0.0056 


0.0053 


XI z 


0.8644 


0.0646 


0.0609 


X2 x 


0.0065 


0.1841 


0.1847 


X2 z 


0.4311 


0.6329 


0.6327 


A3 x 


0.2912 


0.2898 


0.2879 


A3 z 


0.3298 


0.1674 


0.1654 


X4 x 


0.0816 


0.0799 


0.0792 


XA y 


0.5435 


0.5459 


0.5471 


X4 z 


0.1659 


0.3298 


0.3282 



TABLE IV: Structural properties of monoclinic P2\jc 
CsBa2As, A=Br,I. The lattice parameters are from exper- 
iment (Refs. [52] and l53f) . while the internal parameters are 
from total energy minimization using the PBE functional. 





CsBa 2 Br 5 


CsBa 2 I 5 


a, 6, c (A) 


13.816 


9.987 


8.665 


14.637 


10.541 


9.256 
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90.2 






90.194 


Cs 


0.1665 


0.9897 


0.5586 


0.1675 


0.9911 


0.5594 


Bal 


0.1768 


0.5002 


0.5300 


0.1782 


0.4999 


0.5332 


Ba2 


0.4921 


0.2511 


0.5704 


0.4913 


0.2500 


0.5772 


XI 


0.4211 


0.9614 


0.6782 


0.4205 


0.9562 


0.6810 


X2 


0.0980 


0.4602 


0.1643 


0.0956 


0.4565 


0.1669 


A3 


0.4980 


0.2849 


0.9535 


0.4985 


0.2813 


0.9592 


X4 


0.3122 


0.2868 


0.3356 


0.3122 


0.2819 


0.3375 


A5 


0.2192 


0.2732 


0.8099 


0.2216 


0.2694 


0.8110 



compounds occur in a rhombohedral R3 structure, with 
the cations on site 6c (0,0, z) and I on 18/ (u,v,w). With 
the hexagonal setting, the experimental lattice param- 
eters are a=7.4864 A, c=20.880 A, for YI 3 (Ref. \Ej§ 
and a=7.516 A, c=20.7180 A, for Bil 3 (Ref. M). The 
calculated internal coordinates are z=0.1664, w=0.6568, 
v=0.9993, w=0.4156, for YI 3 and z=0.1664, tt=0.0022, 
v=0.6685, w=0.4123, for Bil 3 . 



IV. ELECTRONIC STRUCTURE AND 
OPTICAL PROPERTIES 

As mentioned, three functionals were used in the 
present study, the PBE GGA for the structural prop- 
erties, and the Engel-Vosko and TB-mBJ functionals 
for the optical properties. The recently developed TB- 
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FIG. 1: Calculated band structure of BaF2 with the three 
functionals. 



mBJ has been shown to give much more accurate band 
gaps than other semilocal functionals in a wide variety 
of materials^ As such, we emphasize results obtained 
using that functional. For comparison, the fundamental 
band gaps of the compounds studied are listed in Ta- 
ble [V] The EV gaps are invariably intermediate between 
those of the standard PBE and the TB-mBJ functionals 
with the exception of the La containing compounds. In 
those compounds, the /-resonance is found inside the in- 
sulating gap for all three functionals. The position of the 
/-resonance relative to the valence band edge depends 
on the compound, but for a given material is practically 
unchanged between the different functionals. For these 
compounds we show both the gap between the valence 
band maximum and the bottom of the /-bands making 
up the resonance, as well as the gap between the valence 
bands and the non-/ conduction bands. Experimental 
band gaps are unavailable for most of these materials. 

We start with BaF2 , which as mentioned is a relatively 
well characterized material. The calculated band struc- 
ture with the PBE, EV and TB-mBJ functionals is as 
shown in Fig. [TJ while the corresponding electronic den- 
sity of states (DOS) is given in Fig. H The PBE band 
structure has a direct gap at T that is clearly underesti- 
mated with respect to experiment^ The valence bands, 
of F p character, are narrow with a width of less that 2 
eV, while the metal derived conduction bands are more 
dispersive. The EV band structure is similar to the PBE 
result but with a somewhat larger band gap. The mBJ 
band structure has a considerably enhanced gap of 9.8 
eV. Furthermore, the conduction bands are distorted rel- 
ative to the PBE result, with the conduction band min- 
imum now at the X point, although it should be noted 
that the distinction between direct and indirect is not so 
important for the properties of this compound because 
of the very flat valence bands. 

The experimental band gap of BaF2 is usually quoted 



TABLE V: Calculated fundamental band gaps in eV for the 
various compounds. For La containing compounds we list 
both the band gap from the valence band edge to the La 4/- 
resonance (denoted /) and to the conduction bands excluding 
the resonance (denoted cb, see text). Spin orbit is included 
for Bil 3 . 





PBE 


EV 


TB-mBJ 


BaF 2 


6.87 


7.46 


9.77 


BaCl 2 


5.23 


5.66 


6.45 


BaBr 2 


4.45 


4.85 


5.39 


BalBr 


3.79 


4.13 


4.45 


Bal 2 


3.50 


3.87 


4.14 


Cal 2 


3.49 


3.92 


4.26 


SrBr 2 


4.99 


5.25 


5.91 


YI 3 


2.81 


3.02 


3.31 


Bils 


1.43 


1.67 


1.82 


Cs 2 LiYCl 6 


4.90 


5.30 


5.97 


CsBa 2 Br 5 


4.53 


5.13 


5.77 


CsBa 2 I 5 


3.69 


4.14 


4.49 


K 2 LaCl 5 / 


4.46 


4.46 


4.59 


K 2 LaCl 5 cb 


5.23 


5.49 


6.32 


K 2 LaBr 5 / 


3.69 


3.76 


3.74 


K 2 LaBr5 cb 


4.33 


4.66 


5.26 


K 2 LaI 5 / 


2.79 


2.95 


2.80 


K 2 LaL cb 


3.38* 


3.61 


3.99 


LaCl 3 / 


3.97 


4.04 


4.15 


LaCbj cb 


4.95 


5.26 


5.82 


LaBr 3 / 


3.20 


3.28 


3.24 


LaBr3 cb 


3.98 


4.24 


4.60 



'/-resonance overlaps bottom of conduction band. 



z 




FIG. 2: (color online) Electronic density of states of BaF 2 as 
calculated with the PBE, EV and TB-mBJ functionals. 



as 11 eV based on the measured ultraviolet reflectance 
spectrum^ This spectrum shows several features near 
the gap: a prominent feature at 10.00 eV, assigned as an 
exciton, a weaker feature at 10.57, also assigned as an 
exciton, and a direct band edge at 11.0 eV. Subsequent 
two photon absorption experiments assigned the gap as 
10.6 eV. eo The TB-mBJ direct gap at X of 9.90 eV is 
~ 1 eV lower. Importantly, the experimental spectrum 
has three large peaks in the valence band region below 
16 eV. The calculated TB-mBJ reflectance spectrum is 
shown in Fig. [3] There are two prominent peaks that can 
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FIG. 3: Calculated TB-mBJ reflectivity spectrum of BaF2. 
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FIG. 4: (color online) Calculated refractive index of BaF2 
with the three functionals. 



be matched with the upper two peaks in the experimental 
spectrum (Ref. I59T ) using a band shift of ~ 1 eV. This is 
consistent with the assignment of the lower experimental 
peak as excitonic and the conclusion that the TB-mBJ 
band gap is underestimated by ~ 1 eV. 

The calculated refractive index is shown as a function 




FIG. 5: (color online) Calculated refractive index of 
Cs 2 LiYCl 6 . 



of energy for the three functionals in Fig. [4] The refrac- 
tive indices follow the trend in the band gaps, with the 
result of the TB-mBJ functional lowest. Experimentally, 
the refractive index of BaF2 is practically constant, ris- 
ing from n=1.465 at low energy to n/j=1.474 at 589 nm 
(2.107 eV)i£i This is in close agreement with the results 
using the TB-mBJ functional. We obtain a low frequency 
value of n=1.42 with a weak dispersion as in experiment. 
The low frequency PBE value is higher at n=1.56 and 
shows a stronger energy dependence consistent with the 
smaller gap, while the EV value of n=1.52 is intermedi- 
ate. 

This result shows that the TB-mBJ functional not 
only improves the band gap, but also improves the op- 
tical response of this halide relative to the PBE func- 
tional, which is regarded as state of the art for calcu- 
lations of total energies. This is important because the 
response depends not only on band energies but also on 
the wave functions, specifically through the dipole ma- 
trix elements. The charge density and total energy are 
fundamental quantities within density functional theory 
while the single particle energies are not. It is strongly 
thought that the Kohn-Sham eigenvalues of the exact 
density functional will not produce correct band gaps in 
insulators, even though its charge density and static re- 
sponse will be exacts Therefore it is of interest to note 
that the improvement in the band gap of semi-local func- 
tionals, particularly going to the EV and then the TB- 
mBJ functional improve the static response at the same 
time as they improve the band gap. In the following we 
focus on results obtained with the TB-mBJ functional. 

The calculated TB-mBJ band gap of 5.97 eV for the 
high light output elpasolite scintillator, CS2L1YCI6 agrees 
well with the experimental onset of optical absorption at 
5.9 eV. The refractive index is shown in Fig. [SJ 38 The low 
frequency limit is n=1.65. There is more dispersion than 
in BaF2 , but the spectrum is otherwise featureless almost 
up to the band edge. We are not aware of experimental 
measurements of the optical properties of this material. 

The isostructural tri-iodides, YI3 and B1I3 offer an in- 
teresting contrast. As mentioned, when activated with 
Cc 3+ , YI 3 is a very high light output scintillator, while 
B1I3 is not. The emission of Yl3:Ce 3+ starts at ~ 2.9 eV 
and peaks at 2.26 eV£ The high energy onset is below 
our calculated TB-mBJ gap of 3.31 eV, but is higher than 
the gap obtained with the standard PBE GGA and prior 
calculations done with the local density approximation 
(LDA)^ The difference between YI3 and B1I3 can be 
understood in terms of the electronic structure, specifi- 
cally that the band gap of Bil3 is too small and that this 
is due to low lying Bi 6p states. 

In order for scintillation to take place an electron hole 
pair must recombine at the Ce 3+ site. For this to be ef- 
ficient both the upper and lower states for the activator 
ion should be in the band gap. This is evidently not pos- 
sible in B1I3. As mentioned, the Ce 3+ emission in YI3 
has a high energy onset of 2.9 eV, which is similar to 
that in other iodides. This is a lower bound on the en- 
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FIG. 6: (color online) Comparison of the electronic DOS of 
YI3 and Bi3 as obtained with the TB-mBJ functional. 
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FIG. 7: (color online) Calculated refraction of YI3 (top) and 
Bil3 (bottom) as a function of energy based on the TB-mBJ 
functional. 



ergy difference between the upper and lower states of the 
activator ion. A comparison of the calculated TB-mBJ 
electronic DOS of the two materials is given in Fig. [5] 
The lowest conduction bands in YI3 are formed from Y 
Ad states. In contrast, there is a manifold of Bi 6p states 
making up the conduction bands in B1I3, and these occur 
at lower energy. This leads to a 1.49 eV lower band gap 
in the Bi compound. Our TB-mBJ gap for B1I3 is sub- 
stantially higher than that obtained in prior calculations 



FIG. 8: (color online) Crystal structure of YI3 with the re- 
laxed atomic positions. The large red balls denote I, while 
the smaller gray balls are Y. Note the layered structure. 



using other functionals , 64 ! 65 but it is still only 1.82 eV. 
The experimental optical gap is 2.0 eVj££ This does not 
leave a sufficient energy range for both the upper and 
lower states to be in the gap. This is similar to the prob- 
lem in activating Pb phosphate glasses with rare earths^ 
In fact, considering the fact that, all things being equal, 
lower gaps favor higher light output, but that too low a 
gap prevents activation, the gap of YI3 is probably close 
to optimum for a Ce 3+ activated material. 

The refractive indices of YI3 and B1I3 are presented 
in Fig. [3 Surprisingly, the refractive index of YI3 is 
almost isotropic, even though the material itself is struc- 
turally very anisotropic as shown in Fig. |SJ The low en- 
ergy limits of the refractive indices of YI3 are n xx =2.08 
and n zz =2.06, which is an anisotropy of 1%. B1I3 is 
slightly less isotropic, but is still remarkably so for such 
an anisotropic layered crystal structure. The low energy 
refractive indices 2.76 and n zz =2.7l for B1I3. 

Jellison and co-workers^ reported a value of 3.1 at 1.6 
eV. For comparison, our calculated values at 1.6 eV are 
3.16 in-plane (xx) and 3.02 out-of-plane (zz), in excellent 
agreement with the experiment. 

We find a similar remarkable near isotropy of the op- 
tical properties in several other halides as well. This is 
the case even for monoclinic CsBa2Brs and CsBa2l5 as 
shown in Fig. [9] The three diagonal components of the 
refractive index are practically the same almost up to the 
band edge, while the off-diagonal component associated 
with the monoclinic symmetry is practically zero over 
this energy range. 

For comparison, we also performed calculations for 
monoclinic CdW04. We followed the same procedure re- 
laxing the internal coordinates using the PBE functional 
and then performing optical calculations using the TB- 
mBJ functional. The calculated band gaps are 2.99 eV 
with the PBE functional and 4.16 eV with the TB-mBJ 
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FIG. 9: (color online) Calculated refraction for CsBa2Brs 
(top) and CsBa2l5 (bottom) based on the TB-mBJ functional. 
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FIG. 10: (color online) Calculated refraction for CdW04 with 
the TB-mBJ functional. 



form. The PBE value is similar to that obtained previ- 
ously by Abraham and co-workers*^ The TB-mBJ value 
is close to the experimental value of 3.8 - 4.1 eV (see Ref. 
68h . The refractive index is shown in Fig. [TU] As may 
be seen, in contrast to the halides, it is quite anisotropic 
with a substantial off-diagonal component. 

Returning to halides, we next discuss La containing 
materials. The band structure of these materials may be 
described as that of an equivalent material based on a 




FIG. 11: (color online) Comparison of the DOS with the three 
functionals for LaBr3 (top) and K^LaBrs (bottom). Note that 
the position of the /-resonance that sets the fundamental gap 
is relatively unaffected by the changes in functional. 



simple trivalent ion, e.g. Y, plus additional unoccupied 
/ derived bands that occur in the band gap. The / bands 
comprise the so-called /-resonance. In Table |V] we give 
both the value of the fundamental gap, which is from the 
top of the halogen p derived valence bands to the bottom 
of the /-resonance as well as a larger gap, denoted "cb" , 
which is from the top of the valence bands to the bottom 
of the conduction bands, excluding the /-resonance. 

We find that in contrast to the band gaps of the other 
materials studied, the fundamental gap of these materials 
is practically unchanged upon going from the PBE to the 
TB-mBJ functional. That is the position of the / bands 
relative to the valence band maximum is almost the same 
with these three functionals. This is illustrated for LaBr3 
and K^LaBrs in Fig. [TT] The position of the higher ly- 
ing non-/ conduction bands is, however, increased with 
the TB-mBJ functional. This does lower the optical re- 
fractive index, even though the fundamental band gap is 
unchanged. 

The calculated refractive indices of LaCl3 and LaBi-3 
are shown in Fig. 1121 The low energy values circ ft xx — 

1.87 

and n zz =1.91 for LaCi3 and n xx —2.QA and n zz =2.09 for 
LaBr3. This follows the expected trend where bromides 
have higher refractive index than chlorides. Again, we 
find low anisotropy in materials that are structurally 
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FIG. 12: (color online) Calculated refraction for LaCls(top) 
and LaBr3 (bottom) based on the TB-mBJ functional. 



anisotropic. The same trend with halogen atomic number 
and also very low optical anisotropy is found in K^LaCls, 
K2LaBr5 and K^Lal^, as shown in Fig. [13] The calcu- 
lated TB-mBJ band gap to the non-/ upper conduction 
band edges ("cb") for K^LaChj, K^LaBrs and K^Lals are 
6.32 eV, 5.26 eV and 3.99 eV, respectively, as compared 
to the experimental estimates of 6.6 eV, 5.5 eV and 4.5 

c y 45,46 

The refractive indices of BaCl2, BaBr2, BalBr and 
Bal2 are given in Fig. [TU while those of Cal2 and SrBr2 
are given in Figs. [T5]and[T6l respectively. The calculated 
TB-mBJ band gap for BaCl2 of 6.45 eV is in reasonable 
accord with the estimate from CI- if x-ray spectroscopy 
of 7.0 eV£2- As mentioned, BaCl 2 , BaBr 2 , BalBr and 
Bal2 occur in a complex orthorhombic Prima structure. 
SrBr2 is tetragonal, while Cal2 is rhombohedral. None 
of these materials is cubic. Cal2 activated with Eu 2+ has 
been known to be an extremely high light output material 
since the 1960's.— However, this hexagonal material has 
not been used in applications because of crystal growth 
problems due to its platelet growth habit. Nonetheless, 
again we find only very small anisotropy in the optical 
properties for all of these materials, similar to what we 
found previously for Srl2^ 
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FIG. 13: (color online) Calculated refraction for 
K2LaCl5 (top), K^LaBrs (middle) and K^Lals (bottom) 
based on the TB-mBJ functional. 



V. DISCUSSION 

We have two main conclusions, besides the numeri- 
cal data, which we hope will be useful in improving the 
design of scintillator systems for better light coupling. 
First of all, we find that the newly developed TB-mBJ 
functional greatly improves both the band gaps and the 
optical properties in a broad class of halide materials, 
consistent with results reported for other compounds. 19 
Considering the computational efficiency of this method, 
which is similar to standard density functional methods, 
we expect that this method will enable optical character- 
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FIG. 14: (color online) Calculated refraction for BaCi2, 
BaBr2, BalBr and Bai2, from top to bottom using the TB- 
mBJ functional. 
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FIG. 15: (color online) Calculated refraction for Cal2 based 
on the TB-mBJ functional. 




FIG. 16: (color online) Calculated refraction for SrBr2 based 
on the TB-mBJ functional. 



ization of new complex halide scintillators and perhaps, 
considering that band gap is a key parameter, more ef- 
fective theoretical screens for new scintillators. 

Secondly, and quite unexpectedly, we find that a wide 
variety of halide scintillators based on CI, Br and I are 
practically isotropic from an optical point of view, even 
though many of them are highly anisotropic from the 
point of view of structure and other properties. The 
broad range of materials in which this occurs implies that 
it is a general feature of halide chemistry rather than a 
special coincidence for certain compounds. 

Qualitatively, it may be understood from the coordi- 
nation environments and general band structure features. 
In particular, these materials have relatively wide band 
gaps due to the large electronegativity difference between 
the cations and the halogen atoms. The halogen p derived 
valence bands are narrow compared to typical metal ox- 
ides. This narrow band width suggests that one can un- 
derstand the properties in real space instead of depending 
essentially on detailed band dispersions. Additionally, 
the valence band formation comes at least largely from 
direct hoping between the halogen p orbitals on adjacent 
sites these materials. The structures can be described in 
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anion contact terms,-- From a structural point of view, 
the anion lattices of these compounds are distortions of 
high symmetry structures, with most of the anisotropy 
coming from the cation placement in the intertices of the 
anion lattice. This places the cations in locally highly 
symmetric cages based on high nearest neighbor anion 
coordination numbers. Since the anion bands are rela- 
tively non-dispersive, the main crystal structure depen- 
dence comes from the conduction bands. Therefore, the 
small optical anisotropy of these materials can be ratio- 
nalized in terms of the highly symmetric local environ- 
ments of the cations as far as nearest neighbor anion co- 
ordination is concerned. Further work to elucidate this 
is clearly needed. 

In any case, this result has important implications for 
the development of halide scintillators. Gamma spec- 
troscopy requires high quality uniform crystals of suffi- 
cient size to effectively stop the Gamma rays in the scin- 
tillator volume. This typically requires cm sized crystals, 
and many applications benefit from still larger sizes. As a 
result, crystal growth is one of the main challenges in the 



development of new halide scintillators. The near opti- 
cal isotropy of these materials, however, suggests that 
ceramic scintillators of sufficient size to be useful for 
Gamma spectroscopy can be made. While this will re- 
quire solution of a number of problems, for example, the 
development of methods to produce dense sintered bod- 
ies without contamination in these often hygroscopic, air 
sensitive materials, it may enable the use of low symme- 
try difficult to grow halide materials for gamma spec- 
troscopy and other scintillator applications. Further- 
more, ceramic materials are generally lower cost than 
single crystals, especially if large volumes of material are 
needed. 
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